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As part of a spectroscopic survey of supernova remnant candidates in M83 
using the Gemini-South telescope and CMOS, we have discovered one object 
whose spectrum shows very broad lines at Ha, [O I] AA6300, 6363, and [O III] 
AA4959, 5007, similar to those from other objects classihed as ‘late time super¬ 
novae.’ Although six historical supernovae have been observed in M83 since 
1923, none were seen at the location of this object. Hubble Space TelescopeWide 
Field Camera 3 images show a nearly unresolved emission source, while Chandra 
and ATCA data reveal a bright X-ray source and nonthermal radio source at the 
position. Objects in other galaxies showing similar spectra are only decades post¬ 
supernova, which raises the possibility that the supernova that created this object 
occurred during the last century but was missed. Using photometry of nearby 
stars from the HST data, we suggest the precursor was at least 17 Mq, and the 
presence of broad Ha in the spectrum makes a type H supernova likely. The 
supernova must predate the 1983 VLA radio detection of the object. We suggest 
examination of archival images of M83 to search for evidence of the supernova 
event that gave rise to this object, and thus provide a precise age. 

Subject Headings: galaxies: individual (M83) - galaxies: ISM - ISM: supernova rem¬ 
nants - supernovae: general 

Facilities: Gemini-South Observatory (CMOS), Hubble Space Telescope (WFG3), Ghan- 
dra X-ray Observatory 


Introduction 


The transition from fading supernova (SN) to young supernova re mnant (SNR) is not 
well characterized either observationally or theoretically flFesen l200ll) . For a handful of 
objects where the SN was observed and there have been observations after the SN light has 
faded, so-called ‘late time’ observati ons of the SN positions reveal v ery broad emission lines 
due to the rapidly expanding ejecta flMilisavlievic et al.ll2009l. l2012l. and references therein). 
The visibility of these broad lines at late times seems to be restricted to core collapse SNe 
of various types. Even then, only a small subset of these SNe are detected at late times, 
depending on the mass loss history of the progenitor and/or the conditions in the surrounding 
interstellar or circumstellar medium. 


By way of comparison, spectra of young SNRs for which the SN was not observed can 
show essentially identical spectra, with very broad lines due to rapidly expanding ejecta for 
up to a few thousand years after the explosion. Examples of these ejecta-dominated SNRs 
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include the extraordinary young SNR 

n NGC 4449 ( 

Milisavlievic & Fesen 2008). thought 

to be of order 70 vears post-explosion ( 

Bietenholz et al. 

2010 

), Cassiopeia A in our Galaxy 

fsee Fig. 14 of Milisavljevic et al. (2012 

) for an ‘integrated’ spectrum of this object), which 


dates to ~1680 AD, and even 1E0102-7219 in the Small Ma gellanic Clond, which i s still 
ejecta-dominated eve n ~2000 years post-explosion (IBlair et al.ll2000l: iFinkelstein et ahll2006 


Vogt fc Dopitall201Cl[l . 


Becanse there are relatively few objects available to stndy this evolntionary phase and 
their characteristics are varied, we choose to combine these two gronps together, and consider 
them all to be yonng SNRs. This makes sense from the perspective that there may be SNe 
that are less than 100 years post-explosion bnt for which the SN event was not observed 
(like the NGC 4449 object mentioned above), which makes the distinction rather artificial. 
However, even combining these gronps, there are relatively few objects available to stndy. 
Thns, finding additional members of the class is one way of improving onr nnderstanding of 
the parameter space controlling this transition from fading SN to yonng SNR. 


We are pursning a mnlti-wavelength observational campaign o n the sonth e rn fac e-on 
spiral galaxy M83 (NGC 5236; SAB(s)c; i = 24°; d = 4.61 Mpc; Saha et ah ( 2006 1) to 
stndy a broad range of topic s inclnding especially the young SNR population and the associ¬ 


ated star formation activity flLong et al.ll2012l:lBlair et al.ll2012l:lLong et al.ll2014J:lBlair et ah 


1985 


2014J1. M83 has host ed six recorded supernovae (SNe) since 1923 fICowan fc Branch 
Stockdale et al.l 120061 ). second in number only to another face-on spiral, NGC 6946, with 
nined] The positions of these SNe across the face of M83 are shown for reference in Fig. [T] 
along with the position of the object that is the subject of this paper. 

Optically, SNRs have most often been identified through higher ratios of [S II] AA6716, 6731 
with resp ect to Her than are typically found in photoionize d gas, a signatu r e of sh ock-heated 
emission (IMathewson fc Clarkelll972l) . In the case of M83, iBlair fc Long (12004 henceforth 
BL04) detected 71 SNR candidates using interference hlter imagery in these two emission 
bands and confirmed 20 of these objects with spectroscopy. More recently, we have expanded 
the SNR candidate list to at lea st 225 objects usi ng imagery from the Magellan I 6.5 m tele¬ 
scope and IMACS instrument ( Blair et ah 2012 . henceforth B12)@. Furthermore, we have 
used Hubble Space Telescope Wide Field Came ra 3 (WFC3) images to investigate 63 of the 
smallest diameter candidates (IBlair et al.ll2014fl . about half of which were newly detected in 


^Interestingly, three of the nine SNe in NGC 6946 have o ccurred just s i nce th e tu rn of the millennium : 
SN2002hh, SN 2004et, and SN 2008S. See the SN catalogs by Barbon et al. (1999) and Lennarz et al. ( 2012 1 
for a full accounting in both galaxies. 


^See also IBlair et al.l (l2013l) for a correction to the flux scale reported in B12. 
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the HST survey and many of which are also strong soft X-ray sources flT^ong et ahl I2014J) . 
Between the HST and Magellan surveys, there are now more than 250 SNR candidates in 
M83. 

We are carrying out a new spectroscopic survey using Gemini-South and the CMOS 
instrument to follow-up many of th ese SNR candidate s and elucidate their properties (Win¬ 
kler et ah, in preparation; see also iLong et ahl (12012[) who observed the young remnant of 
SN1957D in M83). As part of this survey, we have obtained the spectrum of an object that 
shows surprisingly broad lines, especially at Ha. The object corresponds to BL04-50 in our 
initial survey and B12-174 in our more recent study, but the object of interest here appears 
as a bright knot or condensation on the northern limb of the larger object. Hereafter, we 
will refer to this bright knot as B12-174a. In this paper, we report an analysis of this object 
using a Gemini South spectrum along with the data from our HST, Chandra, and ATGA 
radio observations, and conclude that it likely represents the remains of a SN that may have 
occurred within the last century but was not reported. 


2. New Spectrum and Other Supporting Data 


In Fig. [2] we show a six-panel context image of the region surrounding B12-174. The 
optical data are from the Magellan imaging program described by B12, where the emission¬ 
line images have been continuum subtracted. I and B continuum bands are also shown. 


along with an aligned three-color image from the Chandra data reported by iLong et ah 


(12014 object X316). The green circle in the image is centered at the position given by B12, 
where the assessment of the object indicated that the bright knot toward the north was likely 
associated with a fainter extended shell of emission to the south and east that not only seemed 
to have elevated [S H]/Ha emission but also emission in [O HI] A5007. Thus, the diameter 
of t he overall object w as measured to be 3" (~66 pc), and the object was not included in 
the iBlair et ahl (12014 assessment of the smallest diameter SNRs. However, in Fig. |2] it is 
clear that the X-ray source is centered on the bright, unresolved optical knot that we now 
refer to as B12-174a, and it is this knot that we observed spectroscopically. Interestingly, the 
fainter extended emission identified as B12-174 seems to also show elevated [S H]/Ha along 
with detectable [O HI] emission, making it also likely to be shock excited. It is tempting to 
suspect an association betw een these two nebulae, as is the case for the LBV nebula B416 
in M33 (IFabrika et ahl 120051) . but a direct relation remains uncertain in the present case. 
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2.1. Spectroscopic Observations 


As part of a Gemini observing program GS-2011A-C-1, Winkler, PI) to obtain 
spectra of SNRs in M83, we obtained a spectrum of B12-174a using the Gemini Multi- 
Object Spectrograph (GMOS) on the 8.2m Gemini South telescope on 2011 April 9 (UT). 
Our observations consisted of three 2000 s exposures, using the 600 line blue grating (G5323) 
and GG455 cutoff hlter, in multi-slit mode with a custom mask. For B12-174a, the slit 
was 1725 wide and 6" long, oriented E-W as shown in Fig. [2l and covered the wavelength 
range 4500 - 7000 A. Reduction proceeded using standard Gemini techniques within the 
IRAhi environment, and flux calibration was achieved using a number of spectrophotometric 
standards observed during the same 3-night classical observing run. We subtracted the sky- 
in the two-dimensional spectrum based on regions with relatively low emission at both ends 
of the slitlet shown in Fig. [2l and then summed over the central 172 to obtain the one¬ 
dimensional spectrum shown in the upper portion of Fig. [3l 


The spectrum shows a blue continuum with several very broad emission features su¬ 
perimposed. As we will discuss below, the continuum is likely associated with blue stars 
within the aperture and is not directly related to the emission line source. Hence, in the 
bottom portion of Fig. |3]we show a version of the spectrum where the continuum has been 
subtracted, leaving only the broad emission lines. The dominant emission line covers the 
entire 6500 - 6800 A range, while fainter broad lines are seen centered at roughly 6300 A and 
5000 A. From this spectrum, it is clear that the object was selected as having an elevated 
[S II]/Ha ratio in the narrow band emission line imagery because both the narrow Ha and 
[S H] hlter bandpasses are contained within the broad 6650 A feature. At 5000 A, the 
emission line is fainter, but because the entire narrow [O HI] imaging hlter passband was 
hlled with emission, the source nonetheless stands out well in continuum-subtracted [O HI] 
imaging data. 


The ~6650 A feature may be dominated by Ha, although it is broader than the other 
lines and shows an extension on the long-wavelength side that could be due to [S II] 
AA6716,6731. Also, give n the strength of [N II] AA6548,6583 with respect to Ha observed 
for other M83 SNRs fsee iBlair et al.l fj2014j) ). the possible contribution of these lines cannot 
be ignored. The structure at the rest position of Ha is a subtraction residual from an over- 
lying narrow Ha component, but much of the structure in the broad line is likely to be real, 
indicating either a kinematically complex emitting structure or possibly an emitting clumpy 


^IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Asso¬ 
ciation of Universities for Research in Astronomy (AURA) under cooperative agreement with the National 
Science Foundation. 
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shell fsee iDanzigerl fl2005[) ). 


The broad ~6300 A line is a blend of [O I] AA6300, 6363 and should be dominated by 
the shorter wavelength line, which is normally three times stronger. However, the line prohle 
shows a nearly flat-topped structure with a narrower peak toward the long wavelength side, 
presumably indicating velocity substructure in the emitting material. 0 The strongest (short 
wavelength) portion of the broad 6650 A feature shows a similar structure to [O I], which 
may provide some indirect evidence that the long wavelength extended structure in this line 
is due to a separate feature, presumably broad [S II] AA6716,6731 (see below). 


The broad line at 5000 A is likely due to [O III] AA4959,5007, which again is normally 
seen in a 3:1 ratio but with the longer wavelength line dominating. This broad [O III] feature 
is faint enough that any expected substructure is lost in the noise, but we can say that the 
feature looks less flat-topped than the other lines ([O I] in particular) and no narrow peak 
is seen on the long wavelength side, indicating possible differences in the kinematics and 
sub-structure of the [O III]-emitting material compared with the other lines. No obvious 
broad H/1 is visible at the current detection level. 


Because of the uncertainties caused by line blending and unknown relative line intensi¬ 
ties, we have chosen a simplihed approach to constraining the intrinsic line widths that may 
be present. We have used the specfit task in IRAF to fit a single-component Gaussian 
for each of the ten lines indicated in Fig. |31 The specfit task varies allo wed parameters 
systematically and uses an associated error array to minimize fjKrisslll994l) . We have hxed 
the central wavelengths of all lines relative to Ha and fixed the relative intensities of the 
[N H], [O I], and [O HI] doublets to values dictated by atomic physics, but all other relative 
intensities were allowed to vary. The red curve superimposed on the lower trace in Fig. [3] 
represents a nominal best fit using a single-component Gaussian for all ten lines, where all 
lines have a width of 5200 kms“^. While the fit to the line profiles is not unique, this 
experiment provides a reasonable assessment of the broad line widths and lends credibility 
to the idea that the red “wing” on the strong Ha-|- [N H] blend is likely due to [S H] with 
a similar velocity profile rather than an extended red wing of Ha emission. This fit also 
highlights some of the most obvious non-Gaussian structure in the broad lines, including the 
relatively narrow spikes centered at ~6383 A and 6662 A. Because of the possibility of line 
blending, it is not clear whether these two features correlate directly in velocity space or not. 


For more informati on on observed fO I] pro files in late time SNe and young SNR spectra and their 


variation with time, see lMilisavlievic et all (l2ni(T[ l. 
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2.2. Context from High-Resolution Imaging 


We ha v e insp ected the HST WFC3 imaging data for this region originally reported by 
Blair et ahl f 2014il . The object is in our Field 3, and was imaged in a number of emission 


line and continuum hlters, as reported in that paper. Using the astrometric scale derived 
for these images, the position of B12-174a is RA(J2000) = 13:37:06.646; Dec(J2000) = 
—29:53:32.61, with an uncertainty we estimate as < 0708. Fig. 0] shows a small, 3" region 
of the HST imagery centered on B12-174a. Indeed, the source is nearly unresolved in the 
subtracted emission line images, appearing essentially the same size as faint stars in the 
continuum panel of the Figure. (However, see below for a more quantitative assessment.) 
The southern extended emission region B12-174 seen in the Magellan images is too low in 
surface brightness to be detected in the current WFC3 data. 


We have used the imexamine task in IRAF to quantify the character of the emission 
line source vis a vis stars in the HST data. Averaging the results for the three emission line 
bands (taking cuts in both rows and columns of the images), we hnd a FWHM of 2.93 pixels 
for B12-174a, while averaging a handful of stars of comparable brightness in the F547M 
continuum exposure results in a FWHM of 2.49 pixels for the stars. Formally deconvolving 
this in quadrature is consistent with an intrinsic size of 1.18 pixels (FWHM) for B12-174a, 
which translates to 0.0466" or 1.02 pc at the assumed distance. Hence, it is clear the size of 
the object must be quite small. 


Interestingly, in the lower right panel of Fig. 01 a faint red (I-band) source is coincident 
with the emission line source in the other panels. While this could be a red star, it seems prob¬ 


able, from comparison with spectra of similar objects in other galaxies flMilisavlievic et ah 


2012|), that some of this red emission arises from broad [O H] A7325 emission (and possibly 


fainter [Ca H] and [Ar HI] emission seen in some other objects) within the F814W filter, 
which turns on sharply at 7000 A. Without spectral coverage, however, this remains specu¬ 
lative. In Fig. |5]we provide an expanded view of the region showing the same color image 
from Fig. 0] on the left versus the U-band image on the right. The yellow circle indicates 
the position of the emission line source, and coincides with the red source, while the green 
circle, offset one pixel (0.04") to the east, is centered on the blue source in the U-band image. 
Hence, a faint blue star lies very close along the line of sight. We note that the vertical extent 
of this Figure is 1725, the same as the CMOS slitlet; hence, it is highly probable that much 
or all of the blue continuum seen in the spectrum in Fig. 0] arises from the nearby blue stars 
visible in Fig. 0]and Fig. |5l and not from B12-174a. 











2.3. X-ray and Radio Context 


B12-174a is coincident with the X-ray source X316 and the radio source A080 from 
Long et ah ( 2014 ). The positions listed by Long et ah ( 2014 ) for the Chandra and ATCA 
sources are both consistent with our HST position to within 0.2". In that paper we list a 0.35 
- 8 keV count rate of 0.343 ± 0.024 x 10“^ cts s“^ and a model-dependent flux of F(0.35 - 


8 keV) = 1.280 ± 0.083 x 10“° photons cm“^ s“^, with a corresponding luminosity of = 
7.27 ± 0.47 X 10^® ergss“^. The source has hardness ratios of (M — S)/T = —0.50 ± 0.08 
and (H — M)/T = —0.06 ± 0.05, where the bands are defined as S =“soft” = 0.35 - 1.1 keV, 
M = “medium” = 1.1 - 2.4 keV, H = “hard” = 2.4 - 8 keV, and T indicates Total. This 
places the object well within t he regime o f other known SNRs although somewhat harder 
than average (see, for example, (Long et al.l (12014 ) Fig. 5 and 6). The X-ray source was also 
listed in the earlier X-ray paper by ISoria fc Wul (120031) . their source 100. The somewhat 
higher count rate in these earlier Chandra data is offset by the changing response between 
2000 and 2011 such that there is no evidence that the source has varied over the intervening 
11 years. 


A closer inspection of the Chandra data from iLong et al.l (120141) shows a total of 316 
counts (0.35-8.0 keV, 272.8 source c ounts and 43.2 ba ckground counts) over the total 729 


ks data set. We used ACIS extract (Broos et ah 


associated response files as described bv iLong et ah 


2010l) to extract the X-ray spectrum and 


(120141) . We fit the source and background 


spectra simultaneously in XSPEC yersion 12.8.1 using common detector and sky background 
models for both spectra. We ht the unbinned spectra and used the C statistic to ayoid the 
well known bias in when there are only a few counts per spectral bin. Once the best-ht 
parameters had been determined using the C statistic, we computed weighted by the 
model to provide an estimate of the goodness of ht. Given the low number of counts, we 


ht the data with single component models for the source. We adopted the IWilms et ah 


(l2000l) model for interstellar absorption (tbabs in XSPEC) and ht with both a no nthermal 
(powe rlaw) model and a thermal model (APEC; see iFoster et ahl (120121) ) in XSPEC (lArnaud 
im . 


The best ht with the nonthermal model resulted in a C statistic of 744.6 with 1043 
degrees of freedom, and a reduced y^ of 1.05. The best-htted value of the Nh was 7.2 x 
10^° cm“^ and the best htted value of the photon power law index was 2.15. The Galactic 
column along the line of sight to M83 is ~ 4.0 x 10^° cm“^, so this derived Nh value seems 
consistent with a modest additional column intrinsic to M83. 


Left unconstrained, our thermal model hts wanted to drive the htted Nh value to 0.0. 
Hence, we forced a minimum value of Nh of ~ 4.0 x 10^*^ cm“^. The best ht thermal model 
resulted in a G statistic of 767.1 with 1043 degrees of freedom, and a reduced y^ of 1.02. In 
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this model the Nh was indeed the minimum value that we allowed in the £t, and the best 
htted value of the temperature was a surprisingly high kT=4.21 keV. Given the already low 
value of the reduced we did not explore multi-component models as they are not justified 
by the quality of the data. 

Statistically, both the nonthermal and thermal models provide acceptable fits to the 
data. For reference, the nonthermal model fit to the spectrum is shown in Figure [6l The 
data and source plus background model are shown along with the background (sky-|-detector) 
model alone to demonstrate at which energies the source spectrum dominates the background 
spectrum. The data have been binned for display purposes only. The nonthermal model pro¬ 
vides a good fit to the data with no systematic pattern in the residuals. We consider the 
nonthermal model to be more realistic, given both the quality of the fit and the reasonable 
value of the derived photon power law index. The fact that the thermal model prefers an 
unphysical column and the best-fit value of the temperature is rather high for a young SNR 
argues somewhat against this model. However, the quality of the spectrum does not allow 
us to exclude the thermal model or a combination of nonthermal and thermal components. 
While the models are somewhat inconclusive, clearly the spectrum is relatively hard com¬ 
pared to those from most other young SNRs. This might indicate the presence of a pulsar 
wind nebula or it might indicate an unusually high temperature for the object. We note a 
similar conclusion was reached for SN1957D by iLong et al.l (120121) . although an even harder 
spectrum and a much higher column density was inferred in that case. 


The source also correspond s with ATCA s o urce A 080 from ILong et ahl fl2014i). a n d had 
been catalogued previously by iMaddox et ahl (120061) as source 41. iMaddox et ahl (120061) 
report 6 cm (4.873 GHz) and 20 cm (1. 446 Ghz) VLA observat ions from 1998 and reprocess 
earlier VLA observation from 1983-84 (iGowan fc Branchlll985l. 6 cm only) and 1990 (6 cm 
and 20 cm) for comparison. The 20 cm flux in the latter two epochs remains nearly constant 
while the 6 cm flux drops rather dramatically over all t hree epochs, from 0.7 5 mJy (1983-84) 
to 0.42 mJy (1990) to 0.31 mJy (1998). Gomparing the Maddox et ah ( 2006 ) 4.873 GHz flux 
of 0.31 ± 0.05 mJy in 1998 to the ATGA 5 GHz flux of 0.19 ± 0.02 mJy i n 2011 appears 


to ind icate that the earlier trend of declining flux has continued at ~6 cm. IMaddox et ah 


(I 2 OO 6 I) also show the steepening spectral index from —0.28 in 1990 to —0.50 in 1998, although 
with substantial uncertainties. 
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Discussion 


3.1. Broad Lines and Line Identifications 


The broad emission lines in the B12-174a spectrum are reminiscent of a small class 
of ‘late time supernovae’ (or very young SNRs from our preferred perspective) that have 
been observed spec troscopically in other nea rby galaxies. The best s t udied of these may be 
SN1979C in MlOO f Milisavljevic et ah 2009), but Milisavljevic et ah 1 2012 ) present data on 
a number of other core-collapse SNe a few decades post-explosion, including SN1970G in 
MlOl, SN1980K in NGC6946, SN1993J in M81, and others. Of the available comparisons, 
the one whose spectrum most resembles that of B12-174a is SN 1 970G in MlOl. Although 
the 2010 spectrum of this object (Fig. 6 of iMilisavlievic et al.l (120121) 1 is noisy, it shows 
broad [O II]A7325, which may also be similar to B12-174a if some or all of the I-band 
emission noted previously is c onfirm ed to be due to [O II]. A number of the other objects 
shown by IMilisavlievic et al.l 1120121) also have strong [O II] A7325. However, the class is 
small and varied, and observed characteristics likely depend on both age and environment. 
Thus, it is not clear th at this similarity to SN 1970G can be used as an age indicator (see 
Immler fc Kuntz 2005). 


The shape of the broad feature near Ha (see Fig. 13]) is much broader and more asym¬ 
metrical than the two other lines of [O I] and [O HI]. Our preferred explanation based on 
the spectral fits discussed above is to assume that the extension on the long wavelength side 
arises from [S H] instead of Ha. Making this assumption allows a single assumed width to 
apply to all of the lines, rather than assuming a highly asymmetrical Ha line. This red side 
extension is not always seen in the spectra of ot her similar objects, althou gh a similar shape 
is present in SN1980K’s spectrum (see Fig. 8 of Milisavlievic et ah 1 20121) ). 


On the other hand, it is not even certain that the 6650 A feature is just Ha. Super-solar 
abun dances are seen in H H regions across the bright optical disk of M83 llBresolin &: Kennicutt 
2 OO 2 I ). and many o ther M83 SNRs with spectroscopy sh ow [N H] A6583 to be significantly 
stronger than Ha IjBlair fc Long 1120041: iBlair et ahl 120141) . The same spectra show the sum 
of both [S H] lines to often be comparable to or stronger than Ha. Hence, if the red side 
extension is really [S H], it could be an indirect indication that a significant portion of the 
broad feature centered near Ha is attributable to [N H] emission in addition to Ha. The 
spectral fits discussed earlier were consistent with the [N H] lines contributing significantly 
to this feature, although the line blending does not allow hard constraints. 

The presence of [N H] and [S H] might also help explain the absence of obvious broad 
H/5 in the spectrum shown in Fig. [3l If Ha dominates the 6650 A peak, significant extinction 
would be needed to reduce H/5 below detectability. This is marginally at odds with the low 
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column densities indicated by the X-ray fits and the generally blue appearance of many 
of the stars in the region. If a significant portion of the 6650 A feature is attributed to 
[N II], a lower extinction would apply. A higher quality spectrum including a more stringent 
constraint on any broad H/? feature would help elucidate this situation. 


3.2. A Very Young SNR? 


If we assume a diameter of 1 pc (consistent with the deconvolved source size inferred 
for the emission source) and an expansion velocity of 5200 km s“^ with no deceleration, 
we derive an age of about 96 years for the SNR. Since many SNe start with much higher 
velocities than 5200 kms“^, it is possible that significant deceleration has taken place, in 
which case the object could be significantly younger than 100 years. If the object has a 
diameter smaller than we estimate from our deconvolution, then again the object could be 
signihcantly younger than 100 years. Given that many of the other members of the class are 
demonstrably much younger than 100 years, we consider it likely that B12-174a arose from 
a SN that occurred during the last century but which was not observed directly at the time 
of the explosion. The detection of the object in the initial SNR su rvey of BL04, using data 
obtained in 1991 April, and the radio flux measurement at 6 cm bv ICowan fc Branchl fll985l) 
from 1983 indicates that the SN must have occurred prior to 1983. The rapid change in 6 cm 
flux between 1983 and 1990 and a slower decline after that may indicate a rapid change in 
local environment conditions near the SN position, possibly indicating a SN that took place 
not long before 1983. 


From available SN catalog information (iBarbon et ahl Il999l: iLennarz et al.N201211 . the 
brightest observed SNe in M83, presumably caught near peak magnitude, are SN1968L 
(Type II, V = 11.9) and SN1983N (Type Ib, V = 11.7). A wide range of core collapse SN 
light curves show drops from peak by 3 - 4 magnitudes over the hrst 100 days flFilippenko 
19971 ). Assuming a SN occurred near the time the galaxy slid into the solar avoidance zone 
and was unobserved for 3-4 months, such a SN could have dropped to V ~ 16 or fainter 
and easily have been missed, depending on the particular SN type and other circumstances. 

Furthermore, a cursory assessment of the historical record readily shows significant 
incompleteness that gets more signihcan t the earlier one goes. Using the online version of 
the Asiago catalog f Barbon et ahl 1999 1^. making a cut at galaxies with velocities < 700 
kms“^ (or roughly within 10 Mpc), excluding known type la SNe, and binning into 20 year 
bins, one sees the following number of SNe per indicated period: 1901-1920: 3; 1921-1940: 9; 


®See http://cdsarc.u-strasbg.fr/viz-bin/Cat7B/sn, 
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1941-1960: 9; 1961-1980: 14; 1981-2000: 17; and 2001 - 2014: 27. Hence, it seems reasonable 
that a SN in M83, especially during the period prior to ~1980, could have been missed. 

It would be important to pin down an exact time since the explosion for B12-174a if 
possible, and so we encourage a search of archival image data at the B12-174a position (see 
Fig. [1]) for any indication of the SN that created this young remnant. Unfortunately, if 
nothing is found, we will not be able to separate the idea that the SN predates the historical 
record from the possibilities that the SN was significantly sub-luminous compared to normal 
core collapse SNe or that it was simply missed. 


3.3. Clues to the Progenitor 


SN1957D, also in M83, forms an interesting comparison to B12-174a. In the late- 
1980’s a spectrum of thi s young remnant showed broad oxy gen lines, in particul ar at [O III] 
A5007 flLong et al.lll989[) . As reported bv iLong et al.l (1199211 and more recently bv iLong et ah 
(120121) . however, this broad emission has faded rather dramatically on a fairly short timescale. 
The hard X-ray source and radio characteristics of SN1957D make a plerion interpretation 
likely for this young remnant. A major difference between SN1957D and B12-174a is the 
presence of broad Ha (and possibly [N H]) in the spectrum of the latter, which hints strongly 
for a type H SN designation for the B12-174a SN. SN1957D was discovered well past maxi¬ 
mum and the SN type was not determined at the time. 


As described in detail bv iKim et al.l (120121) . careful photometry in multiple continuum 
bands can be used to investigate the extinc tion proper t ies ari d characteristic age or range 
of ages in a particular stellar population. iLong et al.l (1201211 used HST photometry and 
this CMD htting technique to determine the characteristics of the stellar population near 
SN1957D. Many of the stars were consistent with a very young isochrone of 4 Myr, but 
the envelope of points seemed to be bounded by the 10 Myr isochrone, indicative of a main 
sequence turn-off mass of at least 17 Mq (and possibly signihcantly higher) for the associated 
cluster of stars. Hence, there is little doubt SN1957D was a core collapse SN of some type. 

We have carried out a similar photometric analysis to estimate the age and main se- 
qu ence turnof f mass for the population near B12-174a, again using the procedures outlined 
bv iKim et al.l (120121) . The object is located within a spiral arm to the SE of the M83 nucleus 
(see Fig. [T]), and appears to be associated with a well-defined grouping of stars, shown by 
the green box in Fig. [71 The arrow in this Figure indicates the position of B12-174a for 
reference. We do not center the region of interest on the B12-174a position, but rather set 
the region by eye to be the grouping of bright stars that seem to be associated. Taking a 



















13 


larger region would only include more ‘background’ stars, while a significantly tighter region 
would just diminish the number of stars available for photometry. The size of the selected 
region is ~55 x 76 pc, which is reasonable for a young, coeval grouping. 


Aperture photometry was carried out in the HST F336W (U), F438W (B), F547M 
(narrow-V), and F814W (I) broadband imaged for 72 stellar objects in the region outlined 
in Fig. [71 The results are shown in the color-color diagram (left) and CMD (right) in Fig. [HI 
We only plot stars with photometric errors less than 0.25 mag in (V - I) color and (V - I) < 
1.2, as redder stars tend to be heavily extincted and d o not have U and /or B colors that are 
required for the correction procedure (see section 2.4 of lKim et al.l f[2012[) b The lines shown in 
both panels o f Fig, p are Padova stellar isochrones (CMD V2.3; http://stev.oapd.inaf.it/cmd; 
Marigo et ah f 2008[ )) for the HSTWFC3 filters, Z = 1.5 Zq, and various ages (6.0 < log (age) 


< 9.95). We applied a Gala ctic foreground extinction of Ay = 0.19 for M83 (from NED, 


Schlafiv fc Finkbeinerl (l201l[B . The spread of the points in the color-color panel of Fig. [8] is 


largely due to differential internal extinction between individual stars. The arrow in each 
panel shows the direction of motion of points in each diagram with increasing reddening. 


Fig. [9] shows the same two plots after applying the star-by-star extinction correction. 
The technique is to unredden the stars in the color-color panel by moving them backwards 
along the reddening vector (toward upper left) unt il they match the isoch rone-predicted val¬ 
ues of the reddening-free Q parameter (see eqn. 1 in Kim et ah ( 2012 ): also Mihalas Sz Binne'^ 
fll98l[) ). The corrected colors and magnitudes are then used in the CMD to det ermine the 
i ntrin sic characteristics of the stars. Again, details are provided in sec. 2.4 of [Kim et ah 

(l2ni2[) . 


There is some spread in the points around the youngest isochrones in Fig. [H but the 
data are best fit with the 7 ~ 10 Myr curves, consistent with the group being old enough 
that no Ha is seen associated with the stars but still quite young. This is very similar to 
the age derived for the association adjacent to SN1957D. Inspection of the isochrone data 
files indicates a nominal MS turnoff mass of 22 - 17 Mq for 7 Myr and 10 Myr, respectively. 
The brightest stars in the corrected CMD are clearly younger than the 15 Myr isochrone, so 
we consider ~17 Mq to represent a reasonable lower limit to the B12-174a progenitor mass, 
again consistent with a core-collapse SN type. The corrected position of the B12-174a source 
in this Figure is somewhat too red, which is consistent with the earlier suggestion that the 
I-band flux for this source may be contaminated by [O II] A7325 emission. 


^Although we do not actually transform to the Johnson-Kron-Cousins UBVI system, we will refer to 
these filters as U, B, V, and I for convenience of presentation in the text. 
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4. Summary 

Adding Gemini-South spectroscopy to earlier Chandra, HST, radio, and ground-based 
optical data sets, we have discovered that the object B12-174a represents a very young 
remnant of a SN in M83 that likely occurred in the past century but was not reported. The 
spectrum shows very broad lines at the positions of Ha, [O I], and [O III] and looks similar in 
many ways to spectra of a handful of ‘old supernovae’ in other nearby galaxies that date to 
only decades after the explosion (the oldest being SN1970G in MlOO). A line-htting exercise 
indicates an approximate expansion velocity near 5200 km s“^ as well as the presence of 
broad [S II] and possibly [N II] in addition to Ha. The presence of a radio source in data 
from 1983 implies a date prior to this for the SN explosion, and we encourage a search of 
archival images of M83 for evidence of any activity at the position of B12-174a. Photometry 
of stars in the direct vicinity in the HST data indicates a young population and a potentially 
massive precursor, and makes an association with a core collapse SN of type H likely. 

Historically determined observational SN rates in nearby galaxies are of interest for con¬ 
straining other more indirect determinations of SN rates, but are fraught with the problems 
of incompleteness and small number statistics. The identihcation of B12-174a raises the pos¬ 
sibility that careful surveys of nearby galaxies may be able to correct these rates for missed 
SNe and arrive at improved SN rates for galaxies such as M83. 


We thank additional colleagues on the companion multi-wavelength surveys of M83 
{Chandra, ATGA, and HST) for useful discussions. WPB and PFW acknowledge travel 
support from the National Optical Astronomy Observatory, as well as STScI grants under 
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respectively. Support was also provided by the National Aeronautics and Space Administra¬ 
tion through Chandra grant number GOl-12115, issued by the Chandra X-ray Observatory 
Genter, which is operated by the Smithsonian Astrophysical Observatory for and on behalf 
of NASA under contract NAS8-03060. WPB acknowledges support from the Dean of the 
Krieger School of Arts and Sciences and the Genter for Astrophysical Sciences at Johns Hop¬ 
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Fig. 1.— A V-band image of M83 from the Magellan 6.5 m telescope with the positions of 
historical SNe marked as plus signs. The location of B12-174a, the object of interest to this 
paper, is shown by a red circle in the arm to the SE of the starburst nucleus. In this and all 
imaging hgures that follow, north is up and east to the left. 
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Fig. 2.— This six-panel figure shows the region surrounding B 12-174 in v a rious bands. 
Panels a-e show Magellan IMACS imaging data as reported by iBlair et al.l (1201211 . with 
panels a, b, and c showing continuum-subtracted Ha, [S II], and [O III] and panels d and 
e showing I-ban d and B-band. Pa nel f shows a three-color representation of the Chandra 
X-ray data from iLong et ahl (120141 ). with red being soft (0.35-1.1 keV), green being medium 
(1.1-2.4 keV), and blue being hard (2.4 - 8 keV) X-ra y s. Th e green circle is 6" in diameter 
and is centered on the position given bv iBlair et al.l (120121) . which assumed the southern 
arc of emission (labelled 174 in the [S II] panel) was part of the same SNR. However, the 
Chandra emission is clearly centered on the bright knot of emission that is offset toward the 
north, which we call B12-174a. The yellow box indicates the position of the Gemini-GMOS 
1'.'25 X 6" aperture used to obtain the spectrum in Figured 
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Wavelength (A) 

Fig. 3.— The upper spectrum shows the calibrated Gemini-South CMOS spectrum of B12- 
174a, extracted from the bright point-like source in the yellow box shown in Fig. [2l offset 
vertically by 0.7 x 10“^^ ergs cm“^ s“^ for clarity of presentation. Very broad emission 
lines are seen superimposed on a blue continuum that is due to other stars within the 
aperture. Below we show the same spectrum after fitting and subtracting the continuum. 
The horizontal scale is observed wavelength, while line identifications are shown offset to M83 
velocity (516 kms“^). The red curve superimposed on the lower trace represents a nominal 
fit using a single-component Gaussian for all ten indicated lines. See text for details. 
































































































20 


Ha subtracted 

■ 

[S II] subtracted 

■ 

[0 III] subtracted 

« ■ * ^ 

' V • 

1” ■. 
i,y,B ^ 


Fig. 4.— The HST WFC3 data for a 3" region centered on B12-174a. The four panels show 
the continuum-subtracted Ha, [S II], and [O III] data, plus a three-color rendition of the 
HST data, where RGB = I-band, narrow V (y) band, and B-band emission. An unresolved 
point source is seen in the emission line images. The faint southern arc of B12-174 seen in 
Fig. [2] is too low in surface brightness to show up in the HST data at this exposure level. 
The yellow arrow indicates the faint red source at the same position as B12-174a. 
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Fig. 5.— An extreme close-up of the HST continuum data near the B12-174a position. At 
left is the same three-color image from Fig. 01 and at right we show the F336W (U-band) 
image for comparison. The green and yellow circles show the differing positions of the source 
in blue and red colors, implying a blue star lies in projection just to the east of the SNR. 
The centers of the two circles are separated by 1 pixel, or 0.04". 
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Fig. 6.— An X-ray power law spectral fit for X316, the X-ray source that aligns with 
B12-174a. The data are indicated by the crosses in the upper panel and the histogram 
represents the model fit to the source. The source model includes a nonthermal component 
(a power law), a sky background component, and a detector background component. The 
lower, dashed histogram is the model fit to the background spectrum that includes sky and 
detector background. The data have been binned for display purposes only. 
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Fig. 7.— A wider (8"3 x 7") view of the three-color HST continuum image, indicating the 
2"5 by 3"5 (or roughly 55 by 76 pc) region over which CMD htting was performed (green 
box). Excess extinction can be seen around the outside edges of the region, especially in the 
east and south. The color scale is the same as used in previous Figures. The arrow shows 
the position of B12-174a. 
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Fig. 8.— Results of JIST aperture photometry for stars in the region surrounding B12-174a. 
(left) A color-color diagram; (right) a color-magnitude diagram. The curves are isochrones 
from the Padova models as described in the text. Each dot is a star from within the rectangle 
shown in Fig. [7] as observed (no correction for extinction). The arrows indicate the direction 
of motion of points due increasing extinction, assuming a standard R = 3.1 extinction curve. 
The point marked with a red ‘X’ is the source at the position of B12-174a. 
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F547M - F814W F547M - F814W 

Fig. 9.— Same as Fig. IH but after applying correction for extinction to each star, as 
described in the text. Most of the stars he on the isochrones for young stars, with modest 
scatter. Again, the point marked with a red ‘X’ is the corrected position of the source 
that is coincident with B12-174a. It appears somewhat too red, consistent with I-band 
contamination by SNR emission, as described in the text. 










